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We present the rigorous microscopic quantum theory of the 
interaction of ultracold Bose and Fermi gases with the electro- 
magnetic field of vacuum and laser photons. The main atten- 
tion has been paid to the consistent consideration of dynami- 
cal dipole-dipole interactions. The theory developed is shown 
to be consistent with the general principles of the canonical 
quantization of electromagnetic field in a medium. Starting 
from the first principles of QED we have derived the general 
system of Maxwell-Bloch equations for atomic creation and 
annihilation operators and the propagation equation for the 
laser field which can be used for the self-consistent analysis 
of various linear and nonlinear phenomena in atom optics at 
high densities of the atomic system. All known equations 
which are used for the description of the behaviour of an ul- 
tracold atomic ensemble in a radiation field can be obtained 
from our general system of equations in a low-density limit. 

I. INTRODUCTION 



In the last decade rapid progress is achieved in laser 
cooling of neutral atoms. This leads to the vigorous 
development of the new field of atom optics. One can 
distinguish linear and nonlinear effects in atom optics. 
Theoretical analysis of the linear effects in atom optics is 
carried out without consideration of atomic interactions. 
On the other hand, it was shown that the consideration 
of dynamical dipole-dipole interactions in the radiation 
field leads to the nonlinear effects such as a nonlinear 
Bragg diffraction and the propagation of atomic solitons. 
That is why in recent years the problem of the interaction 
of ultracold atoms in the field of external laser radiation 
has attracted a lot of attention. 

The analysis carried out by O.Morice et al. (T) and 
J.Ruostekoski and J.Javanainen was mainly concen- 
trated on the investigation of the properties of the laser 
radiation modified by atomic dipole-dipole interactions. 
It was shown that this can be described by the refractive 
index which is governed by the Clausuius-Mossotti rela- 
tion known from classical optics if we neglect quantum 
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correlations. If the quantum correlations caused by Bose 
or Fermi statistics are taken into account, the formula for 
the refractive index contains additional terms defined by 
a position dependent correlation function . 

The modification of the properties of the laser radia- 
tion should have a back influence on the behavior of ultra- 
cold atomic ensemble (for example, the motion of atomic 
beam). An attempt to consider this back influence was 
undertaken by several authors |2j|4|H]. In order to in- 
corporate two-body interactions G.Lenz, P.Meystre, and 
E.M.Wright proceeded by eliminating the transverse vac- 
uum field and replacing it by the phenomenological con- 
tact potential [Q. Dynamical dipole-dipole interactions 
which play more important role than contact interactions 
were taken into account in another papers 10, |-|. How- 
ever, in papers by W.Zhang and D.Walls |3|| and G.Lenz 
et al. H the averaged polarization of ultracold atomic en- 
semble was computed as a function of the incident laser 
field, whereas it should be a function of the macroscopic 
or the local field which are different from the external 
laser field due to dynamical dipole-dipole interactions. 
Wallis Q paid attention to this fact and used the cor- 
rect form of the equations for the electromagnetic field. 
However, his equation for the matter field is written in 
the lower order of the density in comparison with the 
equation for the electromagnetic field ||. This problem 
was considered also by Y.Castin and K.M0lmer Q| and 
J.Ruostekoski and J.Javanainen 0. But the equations 
which are used in these papers are very complicated be- 
cause they are written down in terms of the local field and 
dipole-dipole interactions are presented explicitly in the 
form of the sum over dipole fields. This makes the anal- 
ysis very difficult and the equations, which govern the 
time evolution of the matter fields, remained unsolved. 

In our recent paper j^] we have developed the rigor- 
ous quantum theory of the interaction of an ultracold 
atomic ensemble with the electromagnetic field of vac- 
uum and laser photons. The main attention has thereby 
been paid to the consistent consideration of dynamical 
dipole-dipole interactions. We have shown that the re- 
tardation effects influence significantly the behaviour of 
atomic ensemble in the radiation field. Making use of 
the Lorentz-Lorenz relation, we have obtained the gen- 
eral system of equations for atomic creation and anni- 
hilation operators and the propagation equation for the 
laser field which can be used for the self-consistent anal- 
ysis of various linear and nonlinear phenomena in atom 
optics at high densities of the atomic system. We have 
shown that all equations which are used up to now for 
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the description of an ultracold atomic ensemble subject 
to the influence of the laser radiation can be obtained 
from our general system of equations by series expansion 
with respect to the density. 

In spite of the fact that most of the results obtained in 
H remain valid in the case of Fermi-Dirac gases, we con- 
sidered explicitly only the case of Bose-Einstein statis- 
tics. In addition, our treatment of the electromagnetic 
field was not consistent with the QED canonical quanti- 
zation procedure. The vacuum fluctuations of the photon 
field were separated from the electromagnetic field opera- 
tor. Such an operation leads to the violation of the QED 
canonical commutation relations. The aim of the present 
paper is to fill these gaps. 



II. THE MODEL AND HAMILTONIAN 



+ J dr^(r,t)ftw a V> 2 (r,i) 



kA 



H A i = - J drP(r,t)E in (r,i) , 
Haf = - I drP(r,i)D mic (r,t) 



Here we have introduced the polarization operator P(r, t) 
of the ultracold ensemble 

P(r,*)=d^(r,i)^2(r,*)+£r.c.=P + (r,t) + p-(r,i) , 

(4) 

and the operator of the microscopic displacement field is 
given by 



We discuss the combined system of a gas of identical 
neutral polarizable atoms interacting with an incident 
monochromatic field E iTl (r, t) with frequency . The 
probe field E, n is a classical external field. No operator 
is related with it. It is fixed and therefore not subject to a 
dynamics. For example, it may be a standing or running 
laser field. The quantized center-of-mass motion of the 
atoms is taken into account. The atoms are represented 
as point-like neutral two- level atoms with states |1) and 
1 2) with "bare" transition frequency u> a . The atoms in- 
teract via the resonant dipole-dipole interaction of the 
induced point dipoles. This interaction is mediated by 
the exchange of photons with annihilation and creation 
operators CkA and referring in lowest order to plane 
waves in vacuum with wave vector k and polarization A. 

To take into account quantum statistical effects for the 
two-level atoms, we turn to quantized matter fields rep- 
resented by two-component vectors 



V»(r,t)=Vi(r,t)|l>+^(r,*)|2> . 
The equal time commutation relations are 



(1) 



<Mr),4-(r') 



^l(r),V ,t(r') 



to^V) = <5 ij( 5(r-r') , 
£,.7 = 1,2, q 



0, 



c,a, 



(2) 



where [A, B] c is the commutator of the operators A and 
B, and [A, B] a is the anticommutator, which corresponds 
to Bose-Einstein or Fermi-Dirac statistics, respectively. 
In terms of the operators ipi and ip2 and in the multipo- 
lar formulation of QED |l(]|| the Hamiltonian operator 
takes then the form: 



H = Ha 

n A = t 



Hf + Hai + Haf 



(3) 



dr$t(r,t) 



2m 



ipj(r,t) 



kA 



2irhuJk 



V 



e\Ck\ exp (ikr) + H.c. (5) 



We assume that all contact interactions can be ne- 
glected. In order to find the conditions under which this 
approximation is valid, it is necessary to estimate the ra- 
tio Ud/Ug, where Ud and U g are the mean energies of 
the dipole-dipole interaction and the ground-state colli- 
sions, respectively. Treating the ground-state collisions 
in terms of s-wave scattering and restricting to dipolc 
optical transitions we get the following inequality [pllst] 



^ » 37.5s , 



(0) 



where s is the saturation parameter of the atomic tran- 
sition. Thus, the ground-state collisions are negligible 
if the saturation parameter s is of the order of 0.01 or 
higher. This shows, that we don't need very high values 
of the parameter s in order to neglect contact interac- 
tions. 



III. HEISENBERG EQUATIONS OF MOTION 
FOR THE ATOMIC AND PHOTONIC 
OPERATORS 



Making use of the Hamiltonian (||) and the commu- 
tation relations (||) we obtain the following Heisenberg 
equations of motion for the atomic and photonic opera- 
tors: 



if) 



dt 



h 2 V 



^i(r,i)-dEr(r,^ 2 (r,f) 



2m 
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kA 



2 



htp 2 (r,t) ^2 9k\ exp (ikr) c k A (*) , (7) 
h 2 V 2 



kA 



2m 
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kA 



(8) 
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dck\(t) 

at 



h$i(r, t) X ^kA exp (ikr) CkA(i) • 

kA 

ftw k c kX (t) - %kA J dr ex P 
$(r,t)&(r,t)+$(r,t)&(r,t)l , (9) 



.9kA 



12-KUJk , 



where E in are the positive and negative frequency parts 
of the incident classical electric field. The operator prod- 
ucts in Eqs.(Q),(||),(^) are taken in normally ordered 
form. 

The formal solution of ([)]) for the photon operators is 
Ckx(t) = CkA(O) exp (-iukt) 

+*5kA J dt' J dr' exp [iwk(t' — t) — ikr'] 
4>i(r\t>)Mr',t') + ti(r\t')Mr',t')} , (10) 

where the first term 6kA(0) refers to the free-space photon 
field and the second one goes back to the interaction with 
the atoms. 

To study the back reaction of the photons on matter 
we insert ( [l0| ) in (R) and (j^). In the rotating wave ap- 
proximation we obtain the dynamical equations for the 
operators of the two matter fields 



ifi 
1T1 



at 

dj>2(r,t) 
dt 



2V72 



rv 

2m 

n 2 v 2 



^(r,t)-dEj^(r,t)^(r,t), (11) 
-4> 2 {v,t) + h (uj ci + S- 1'7/2) ^a(r, t) 



2m 



(12) 



where S and 7 are the Lamb shift and the spontaneous 
emission rate of a single atom in free space, respectively, 
and we have introduced the operator of the local electric 
field 

E+ c (r,*)=E+(r,i) 



kA 



El 2-K%U>k 
\j — — — e A c k A(0) exp (ikr - iw k t) 



V 



and Vx refers to the point r. Note that in eq.(p"3|) a 
small volume around the observation point r is excluded 
from the integration. 



Eq. ( |l3| ) shows that E^ c (r, t) is the superposition of 
the incident field E^(r,i), vacuum fluctuations of the 
photon field, and the electric field radiated by all other 
atoms. As one should expect it is this local field, which 
drives the inner atomic transition in (|Tl|), (|l2|). 

Equations (|ll|), ( |l2|) can be regarded as an atom- 
optical analogue of the optical Bloch equations H,0. 
They describe the dynamical evolution of the second 
quantized matter in the field of electromagnetic radia- 
tion. 



IV. LORENTZ-LORENZ RELATION AND THE 
SYSTEM OF MAXWELL-BLOCH EQUATIONS 
IN ATOM OPTICS 



A. Local-field correction 



Treating the quantum g medium, it is impor- 

tant to establish a connection to the observables of the 
quantum electrodynamics in media. Having an applica- 
tion to an ultracold gas in mind, we have to relate the 
local field E; oc (r,i) to the macroscopic (or mean) field 
E mac (r, t) which is obtained by averaging in space over 
a region which contains a great number of atoms. For 
a very low density ultracold gas the region on which the 
macroscopic field is averaged may contain actually fewer 
than one atom, yet the averaging technique is still valid 
due to atomic derealization. 

We generalize the idea of Lorentz as it is described in 
p] |l3[ and assume the local field to be a sum of several 
contributions: 

E ioc (r, t) = E„ ear (r, t) + E mac (r, t) - E P (r, t) . (14) 

Let us consider a macroscopically small, but microscop- 
ically large volume V around r. E near is a contribution 
from the dipoles within V. E mac is the macroscopic field 
obtained from the field of all dipoles in the medium. Be- 
cause the dipoles in V are already taken into account in 
E ne ar, we have to subtract from E mac the contribution 
Ep of these dipoles as it is obtained in the averaged con- 
tinuum approximation. Ep is therefore a macroscopic 
quantity. It can be related to the polarization, which is 
an averaged quantity too, according to 



Ep(r,t) 



4tt - , , 
• T P(r,t) 



(15) 



For a novel derivation see [ pT[ . In the following we re- 
strict to the purely classical local-field correction which 
corresponds to a vanishing E„ ear |11-13|. Thus, we see 
that, because of the influence of the near dipole-dipolc 
interactions, the local field is obtained from the macro- 
scopic field in adding the local-field correction fD] [13[ . 



3 



E± c (r,t)=E?. 



(16) 



This equation is also often called the Lorentz-Lorenz re- 
lation. The propagation of the operator of the macro- 
scopic field is thereby given by Maxwell wave equation 
for a charge-free and current-free polarization medium. 



B. Nonlinear matter equation 



this we drop in the following the vacuum fluctuations and 
the spontaneous emission rate 7 from our equations. 

Then substituting (|2l|) in ( |l8| ) we obtain as the in- 
tended result the nonlinear dynamical equation for the 
matter field operator ijji(r,t) 



if) 



dt 



h 2 v 2 



2m 4A 2 (r,t) 



ft + (r) 



(22) 



In the next step we decouple the dynamical equations 
for the two matter field operators. We substitute ( |l6| ) 
in (|ll]) and ( |l2| ) and pass to the reference frame rotat- 
ing with the frequency ujl of the incident field which is 
assumed to be monochromatic 



&(r,t) 



to obtain 



ifl 
■fi 



dt 
d(j>2 
dt 



■ipi 



2m 
S 2 V 2 



2m 

-fi(A + i 1 /2)4> 2 , 



4>2 (r, t) exp (-iu>Lt) , 



fi- , - 4i „ - t - - 



2 --^fi+(r)- 



47T 



(17) 



(18) 



(19) 



with the detuning A = a;^ — w a — <5. The position depen- 
dent Rabi frequency ± (r) = 2d£^ lac (r)/h is related to 
the macroscopic electric field. 

We assume sufficiently large detuning such that the 
contribution of the spontaneous emission is small. We 
may therefore apply the adiabatic approximation to (19), 
which gives 



02 M) 



fi+(r)ViM) 



2 A l (r,t)+i 1 /2 



(20) 



and eliminate the field ^2 from our scheme. The position 
dependent local detuning introduced in 



A i (r,<) = A + ^d 2 ^(r,t)V3i(r,t) 



(21) 



is an operator depending on the density operator of the 
ground state. We interpret this as a shift of the inter- 
nal energy levels which may increase or decrease with 
increasing density depending on the sign of the detuning 
A. Analogous frequency shifts can be found in nonlinear 
optics (see, for instance, |lT| , |l5|| and references therein). 

Because we are mainly interested in atom optical prob- 
lems and want to study the coherent evolution of the 
center-of-mass motion of the gas, we shall neglect spon- 
taneous emission. This is valid for situations where the 
absolute value of the local detuning is much bigger than 
the spontaneous emission rate |A;| ^> 7. In order to do 



It couples to the macroscopic electric field via the Rabi 
frequency fi ± (r). 

On the assumptions stated in the Section I, equation 
(^||) is general. No approximation apart from the adia- 
batic approximation and the rotating wave approxima- 
tion has been made. The dynamical dipole-dipole in- 
teraction mediated through the exchange of photons is 
completely contained. The nonlinearity goes back to the 
density dependent local detuning A;(r, t) of eq. (pl|). For 
increasing density and positive detuning A the local de- 
tuning grows and correspondingly the nonlinear term in 
( p2| ) representing the coupling to the macroscopic elec- 
tric field becomes smaller. On the other hand, for neg- 
ative detuning the absolute value of the local detuning 
decreases with the increase of the density and the non- 
linearity becomes greater. 



C. Microscopic and macroscopic Maxwell equations 

In the multipolar formulation of QED the operator 
of the microscopic displacement field is given by eq.(|5|) 
and the operator of the magnetic field has the form 



B mic (r,i) =i^2\/ 2lT ^ k ¥h (k x e A )c kA (t)cxp (ikr) 



kA 

-H.c. 



V c 



(23) 



Calculating time and spatial derivatives of the opera- 
tors (Eh, (H|) we get microscopic Maxwell equations 

1 d - 

V X E m ^ c 777^mic7 divH m ic 0; 

c at 
19 . 

V x H mic = -— D mic , divDmie = 0, (24) 

c at 



where D r 



E, 



47rP and B r 



The 



same equations were obtained by G.Juzeliunas [p.6| in a 
polariton model of a dielectric medium interacting with 
light. 

In order to get the macroscopic Maxwell equations we 
have to average eqs.(|24|) over the microscopically large, 
but macroscopically small volume of space. However, 
in the case of ultracold atomic gases due to the atomic 
derealization any point of space contains contributions 
from all the atoms. In this sense the microscopic field is 
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already averaged and we may simply replace the micro- 
scopic fields by their macroscopic analogues. 

Assuming that the spatial variations of the atomic den- 
sity are not rapid and taking into account the adiabatic 
solution (pQ), one can rewrite the system of macroscopic 
Maxwell equations in the form of the wave equation 



2<?± 



(25) 



where the refractive index is given by the Clausius- 
Mossotti relation 



1 + ^o^i 



1 



(26) 



Eqs.(p2[),(|25|) and ( |26| ) constitute the general system of 
Maxwell-Bloch equations, which can be employed for the 
solution of different problems in nonlinear atom optics. 
In our paper we have shown that the equations which 
are usually used for the description of an ultracold en- 
semble put in an external laser field can be derived from 
( p2| ) , j25| ) and (^6|) assuming that the density of atoms is 
a small parameter. 



In eq.(pO|) the 5- function vanishes, if ri ^ r 2 . If ri = 
r 2 , the term in square brackets in the r.h.s. of eq.(|30|) 
vanishes. Therefore, the commutator ( |30| ) vanishes for 
any ri and r 2 , and we have the canonical commutation 
relations for the operators of the local electromagnetic 
field. 

Starting from the Lorcntz-Lorenz relation ( [i"6| ) and 
making use of the equal-time commutation relations (|3C|), 
( pl| ) we get the canonical commutation relations for the 
operators of the macroscopic electromagnetic field: 



Erj 



3 ( ri ,t),E™ ac (r 2 ,i) =0 



E moc( r ii t), B^ ac (r 2 , t) 



iAirhc 2J £ mn i-Q^-S(R) 



(32) 



(33) 



The commutation relations (p0[)- (p3|) show that our 
theory is consistent with the general principles of the 
canonical quantization of the electromagnetic field in a 
medium. 



D. Equal-time commutation relations for the 
operators of the electromagnetic field 



V. CONCLUSION 



The operators of the local electric and magnetic field 
strength have the following form: 



E; oc (r,i) 



B ioc (r,i) 



kA 

8tt 



El 2-KflUJk 
\j — 77 — exCkxit) exp (ikr) + H.c. 



V 
P(r,i) 



(27) 



'■^2\l—TT^~r ( k x e A) Ck\{t) exp (ikr) 



kA 

+H.c. 



V 



(28) 



Taking into account the commutation relations for the 
photons and matter-field operators 



= 1: 



0. 



(29) 



we obtain the following equal-time commutation relations 
for the operators of the local electromagnetic field 



87T 



d m d n 5(v 1 -v 2 ) 



$ (r i , t)& (r 2 , t) ~ ft (r 2 , t)fa (r x , t) 

(n , t)$ 2 (r a , t) - i>t ( r 2 , t)i> 2 (n , t) 



(30) 



Er o l c (r!,t),Br oc (r 2 ,t)j = -i^?icJ2 e mn i-^-6(K). (31 



We have treated the interaction of dense Bose and 
Fermi gases with light in the microscopic approach. Our 
results are not limited to low densities and include the 
dipole-dipole interaction consistently. The present anal- 
ysis is concentrated on the common features of the in- 
teraction of Bose and Fermi gases with the radiation 
field. The differences appear when we consider correla- 
tion functions which reflect Bose or Fermi statistics. The 
effects caused by the statistics where considered in , 
where it was shown that the refractive index contains ad- 
ditional terms. However, these statistical effects are not 
very important because the leading term is governed by 
the Clausius-Mossotti formula. 
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